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We demonstrate that quark matter in the color-flavor 
locked phase of QCD is rigorously electrically neutral, de- 
spite the unequal quark masses, and even in the presence of 
an electron chemical potential. As long as the strange quark 
mass and the electron chemical potential do not preclude the 
color-flavor locked phase, quark matter is automatically neu- 
tral. No electrons are required and none are admitted. 



If cold, dense quark matter is approximated as nonin- 
teracting up, down and strange quarks, it has long been 
understood that equilibrium with respect to the weak 
interactions together with the relatively large mass of 
the strange quark imply that the strange quarks are less 
abundant than the other quarks. Thus, noninteracting 
quark matter is electrically positive and a nonzero den- 
sity of electrons is required in order to obtain electrically 
neutral bulk matter. 

Explicitly, weak equilibrium imposes 
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where 3/1 is the chemical potential for baryon number 
and /x e is that for electron number. In the absence of 
interactions, the corresponding number densities are 
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(Throughout, we set m u = rrid — m e = and m c = 
nib = m t — oo.) Electric neutrality requires 

jN u - ±N d - ±iV s - N e = . (3) 

With m s — 0, this can be satisfied with N u — Nd = N s 
and N e = 0. Because m s > 0, however, noninteracting 
quark matter must have fi e > and an electron density 
N e > if it is to be electrically neutral. The condition (1) 
cannot be modified by the strong interactions among the 
quarks. It has long been known, however, that interac- 
tions can modify the relations (2) [1]. We argue here that 
interactions modify these relations qualitatively: they fa- 
vor a state of quark matter in which N u — Nd — N s and 
N e = even when m s ^ and fi e ^ 0, as long as neither 
is too large. 

It is becoming widely accepted that at asymptotic den- 
sities the ground state of QCD with m s = is the color- 
flavor locked (CFL) phase [2-4]. In this phase, color 



gauge symmetry is completely broken, as are both chiral 
symmetry and baryon number (i.e., the material is su- 
perfluid). The effective coupling is weak (because QCD 
is asymptotically free), and the ground state and low- 
energy properties can be determined by adapting meth- 
ods used in the theory of superconductivity (BCS theory) 
[2-5]. The CFL phase persists for finite masses, and even 
for unequal masses, so long as the differences are not too 
large [6,7]. It is very likely the ground state for real QCD, 
assumed to be in equilibrium with respect to the weak 
interaction, over a substantial range of densities. 

Since the CFL phase occurs at weak coupling, it might 
seem natural to think that, like noninteracting quark 
matter, CFL quark matter with m s ^ is electrically 
positive, and that its existence in bulk requires a neutral- 
izing electron fluid. Indeed, this has been tacitly assumed 
in the literature. The presence of an electron fluid dras- 
tically affects the low-energy dynamics of dense matter 
in the CFL phase. Specifically, for example, the elec- 
tron fluid dominates the low-temperature specific heat 
and powerfully resists the motion of magnetic field lines. 

In this note we demonstrate that in reality the CFL 
phase requires equal numbers of u, d, s quarks, and is 
therefore automatically electrically neutral. No electrons 
are required. None are present, even when /i e is nonzero. 

There is precedent for such behavior. In an ordinary 
superconductor, one may consider the effect of a pertur- 
bation SC — (5/t(e|e| — ejejj that splits up and down spin 
energies. It has long been known that for small 5/j., the 
superconducting ground state is completely unchanged 
by such a perturbation and, in particular, it contains 
equal numbers of up and down spins [8] . The analogous 
phenomenon was recently analyzed in QCD, in the con- 
text of pairing between two flavors of quarks with chem- 
ical potentials 



/Hi = fj,- <fyt; /i2=M + ^M 
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For < 5n < with Ao the superconducting 

gap, the ground state of the system is precisely that ob- 
tained for <5/i = [9,10]. By introducing m s ^ in the 
Lagrangian, we are considering a rather different type 
of perturbation, varying the relative mass of the paired 
components, which changes (for example) the relative ve- 
locities within a pair. This perturbation does change the 
form of the superconducting ground state; nevertheless, 
the number of particles of different types remains equal. 

The argument relies on the fact that the CFL ordering 
involves pairing between quarks with equal and opposite 
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momenta and with different flavors. For such pairing 
to be maximally effective, the Fermi surfaces for differ- 
ent flavors of quarks must occur at equal magnitudes of 
the momentum (as opposed to energy). To put it vividly, 
the Fermi surfaces are rigidly locked in momentum space. 
Since the number of occupied states, given by the occu- 
pied volume in momentum space, is now the same for 
all quark flavors, so is their number density. Deviation 
from N u = Nd = N s would reduce the free energy in the 
absence of interactions; CFL pairing, however, reduces 
the free energy most strongly if N u = Nd = N s , and this 
equality is therefore enforced. 

We now make the argument concrete. We work in a 
model in which quarks interact via a four-fermion in- 
teraction which we take to be that with the quantum 
numbers of singlc-gluon exchange [4]. The argument is 
sufficiently general, however, that it applies qualitatively 
(and quantitatively if A /p, is small) to any model with 
four-fermion interactions which are attractive in the ap- 
propriate channel and to QCD at asymptotically high 
density, where quarks interact by single-gluon exchange. 

We can demonstrate the physics of interest by fo- 
cussing on pairing between, say, red up quarks and green 
strange quarks. Let us call these two species of quarks 
"1" and "2", assume they have masses mi = and 
m 2 = m s , and denote their chemical potentials as in (4). 
The generalization of the derivation of the gap equation 
for Ao given in Section 4.3 of Ref. [4] to the case with 
m s 7^ and 5/j, ^ is straightforward, although the al- 
gebra is somewhat involved. The result is 



Ao 

4G 



f d*p 
J (2^) 4 



A w 



(4/x 2 - m 2 s )p 2 - (m + ipo 



Sfi) 2 m 2 ' 
(5) 



where w = Aq+^ 2 +p 2 + (po+i6fi) 2 , wherep= \p\, where 
the four-fermion coupling constant G is normalized as in 
Ref. [4] , and where we have chosen to work in Euclidean 
space. This gap equation has been derived and solved 
explicitly with S/j, = and m s ^ in Refs. [6,11] and 
with m s — and 5fi ^ in Ref. [10]. All we will need, 
however, is the fact that upon replacing the integration 
variable po by p' Q = po + iSfx, the gap equation is seen 
to be independent of #m- This argument holds as long 
as the shift from po to p' does not move a pole from the 
upper-half plane to the lower-half plane, or vice versa. 

We wish to obtain the number densities by differenti- 
ating the free energy density f2 with respect to /x, and so 
must construct f2. For nonintcracting quarks with Fermi 
momenta v\ and v 2 , 

1 [ Vl 

^frcc(Ml,M2,^l,^2) =— P 2 {p-^l)dp 
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+ ^2 [J P 2 (VP^ + ^f - M2) dp. (6) 
We assert that the free energy of the BCS state is 



^BCS = ^frcc (Ml)M2) V\, V 2 ) 
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where the "integrand" is that on the right hand side of 
the gap equation (5), where A solves (5), and where 
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That is, we assert that the correct way to construct the 
BCS state is to first fill noninteracting quark states up 
to the common Fermi momentum (8), and to then pair. 
(The last term in (7) is the condensation energy [4].) 
Why is the procedure embodied in (7) correct? Thinking 
of v\ and v 2 as variational parameters, under what cir- 
cumstances does (8) minimize the free energy (7)? Note 
first that if v\ = v 2 = v is imposed, then the choice (8) 
for v minimizes f2f rco , and therefore minimizes SIbcs be- 
cause the condensation energy is independent of v. Now, 
what about variation with respect to 5v = v 2 — v\1 Try- 
ing b~v > reduces f2f roe , but exacts a cost in reduced 
condensation energy: if v 2 > v\ there is a region of mo- 
mentum space v\ < p < v 2 wherein pairing is impossi- 
ble. (The pair creation operator tries to create either one 
quark of each type or one hole of each type, and there 
is already a "l"-quark and a "2"-hole at every point in 
this region of momentum space. See Ref. [10] for further 
analysis of such "blocking regions" . ) To lowest order in 
m 2 /p, 2 , Sfi/p, and Ao/M, the free energy cost of increasing 
8v from zero outweighs the free energy benefit if 
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This condition can alternatively be derived via analysis 
of the locations of the poles in the integrand in (5): the 
requirement (9) is the requirement that shifting po — > p' a 
in (5) not move a pole across the real axis. 

The condition (9) has a simple interpretation: when 
the free energy gained either by converting a strange 
quark near the common Fermi surface into a light quark 
(m 2 /2jl — 2<5/i) or by converting a light quark into a 
strange quark (25fi — m 2 /2p) compensates for the free 
energy lost by breaking a single pair (2A ), the paired 
state is unstable. Eq. (9) is the criterion for the exis- 
tence of the BCS phase as a local minimum of the free 
energy. To check that it is the global minimum, we must 
compare Qbcs to that for the unpaired state: 
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We use the gap equation to eliminate G in (7), work to 
lowest order in m 2 /p, 2 , S/j/jl and Aq/m, and find 



f^BCS — ^normal 
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meaning that the BCS state is the global minimum of the 
free energy if 



4/z 
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(12) 



Wherever (12) is an equality, there is a first order phase 
transition from the BCS state to the unpaired state. At 
this transition, the Fermi surfaces relax to the (sepa- 
rated) values favored in the absence of interaction. For 
m s = 0, this agrees with previous results [10]. 

Wherever (12) is satisfied, and the paired state is fa- 
vored, we now use the fact that the gap equation is inde- 
pendent of S/j, to conclude that SIbcs of (7) is independent 
of 5fi and therefore 
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To lowest order in A //2, 
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Thus, N\ = N 2 in the paired state even when m s =/= 
and S [i 7^ 0! Note that A , N\ and N 2 all depend on 
m s . The point is that Ni — N 2 , independent of m s . The 
dependence of Ao on m s has been analyzed previously 
[6,7]; the only reason that these authors failed to notice 
that Ni = N 2 is that they did not calculate Ni and N 2 . 

The complete analysis of the CFL state with m s 7^ 
requires the 9x9 block-diagonal color-flavor matrices 
given in Ref. [6]. The analysis is more involved, but 
the conclusion generalizes: the number densities of all 
nine quarks (three colors and three flavors) are the same 
in the CFL phase, even when m s and 5/j, — [i e /2 are 
both nonzero. The excitations in the CFL phase include 
charged Nambu-Goldstone bosons [5], but this does not 
change the analysis of electrically neutral bulk matter: 
adding equal numbers of electrons and positively charged 
mesons costs free energy and is not favored. We conclude 
that quark matter in the CFL phase is electrically neu- 
tral in the absence of any electrons. Even an imposed /x e 
cannot push electrons into the quark matter, because in- 
troducing electrons while maintaining charge neutrality 
and weak equilibrium costs too much pairing energy. 

As an example, take m s — 200 MeV and consider 
quark matter with fi — 400 MeV and fi e = 25 fi = 
150 MeV. (Contact with ordinary nuclear matter would 
impose \i e > 0.) According to (12), this quark matter, 
even with such a large /x e , is in the electron- free CFL 
phase as long as (A„/V2) > [25 - 75 1 = 50 MeV. The 
value of A is uncertain, but A ~ 100 MeV is not un- 
reasonable [4] . Because the stresses imposed on the CFL 
phase by m s 7^ and by [i e > have opposite sign, it 
is more likely than previously thought that, if present, 
quark matter within neutron stars is in the CFL phase. 



The criterion which defines the region wherein the CFL 
phase is favored will deviate somewhat from (12). First, 
there is now an electronic contribution to fi norma i (but 
not to ObcsO- This contribution is only of order 
Second, instead of comparing the CFL phase to a phase 
with no pairing, we should compare it to phases with 
less symmetric pairing. The CFL vs. 2SC comparison of 
Ref. [6] does this at jj, e = 0, and the value of m s at which 
the unlocking transition occurs is in good agreement with 
(12), demonstrating that, at least at /i e = 0, this is not a 
large effect. Third, we expect that, as has been demon- 
strated at m s — [10], there is a region of (m s ,yu e ) just 
outside the CFL region (12) where crystalline color su- 
perconductivity, the analogue of the Larkin-Ovchinnikov- 
Fulde-Ferrell (LOFF) state in ordinary superconductivity 
[12], is favored. Outside the CFL region, very weak pair- 
ing among like flavor quarks is also possible [13]. Because 
both LOFF and single-flavor condensates have much less 
condensation energy than the CFL phase; their effect on 
the location of the boundary (12) is negligible. 

We have demonstrated that while CFL ordering is 
maintained, there will be strictly equal numbers of the 
three types of quarks and rigorous electrical neutrality, 
in the absence of any electrons. If m s or fi e becomes 
too large, however, some less symmetric phase of quark 
matter will have lower free energy. A first order phase 
transition occurs at this boundary as the rigidly locked 
Fermi surfaces spring free under accumulated tension. 

The enforced neutrality of the color-flavor locked phase 
has many consequences: 

In the CFL phase, there is an unbroken U(1)q gauge 
symmetry and a corresponding massless photon given 
by a specific linear combination of the ordinary pho- 
ton and one of the gluons [2,14]. With respect to the 
Q-charge associated with this unbroken U(l) symmetry, 
CFL quark matter is electrically neutral at zero tempera- 
ture, and is a perfect insulator. Because of the absence of 
Q-charged excitations, CFL quark matter is transparent 
to Q-photons and any Q-magnetic flux can move unim- 
peded. In contrast, because electrons have nonzero Q- 
charge, if they were present they would scatter Q-photons 
and would make the material a very good conductor in 
which Q-magnetic flux would be frozen in place [14]. The 
absence of electrons therefore changes the conclusions of 
Ref. [14]: quark matter in the 2SC phase in the core 
of a neutron star anchors magnetic fields, as described 
in Ref. [14]; quark matter in the CFL phase, however, is 
electron-free and therefore offers no resistance to the mo- 
tion of Q-magnetic flux as the neutron star spins down. 

Similarly, if electrons were present they would domi- 
nate the specific heat, which plays a role in the cooling 
of neutron stars [15,16]. In the absence of electrons, the 
specific heat at zero temperature is that of a neutral su- 
pcrfluid, much less than previously thought. The quali- 
tative conclusion that the cooling of a neutron star with a 
CFL quark matter core is dominated by the (large) heat 
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capacity of the nuclear matter mantle remains, however. 

If neutron stars have CFL cores, the absence of elec- 
trons and consequent reduction in specific heat and in- 
crease in transparency amplifies effects (described in 
Ref. [17]) imprinted on the time distribution of the neu- 
trinos from a supernova by a transition from quark-gluon 
plasma to CFL quark matter as the hot, seconds-old pro- 
toneutron star cools. Effects of the first order nature of 
the transition need further investigation, however. 

Schafer realized that if electrons were required to main- 
tain charge neutrality, an alternative would be a conden- 
sate of negatively charged kaons in the CFL phase [18]. 
With N u — Nd = N s in the CFL phase, however, we 
need neither kaon condensation nor a fluid of electrons. 

The broader lesson is that at temperatures which are 
nonzero and small compared to Ao, the transport and re- 
sponse properties of CFL quark matter, in the real world 
with nonzero m s , are dominated by the lightest excita- 
tions of the CFL quark matter itself, and not by electrons 
as had previously been assumed. These bosonic degrees 
of freedom are the massless neutral Nambu-Goldstonc 
boson associated with spontaneous baryon number viola- 
tion (superfluidity) and the neutral and charged pseudo- 
Nambu-Goldstonc bosons associated with spontaneous 
chiral symmetry breaking [2] , which have masses of or- 
der ^/m s m M jA /fi, of order ten MeV [5,4]. The effective 
field theory which describes these light degrees of freedom 
(and thus, we now sec, the phenomenology) is known and 
at high enough density all coefficients in it can be deter- 
mined by controlled, weak-coupling calculations [5,4]. 

Finally, the transition from an ordinary nuclear mat- 
ter mantle to a quark matter core at some radius within 
a neutron star may be greatly simplified if the transi- 
tion occurs directly to quark matter in the CFL phase. 
With a noninteracting quark matter core, one has to face 
the fact that at any given jl, electrically neutral nuclear 
matter and electrically neutral quark matter generically 
have different values of fj, e . Since fj, e must be continu- 
ous across any interface, a mixed phase region is thought 
to form, within which positively charged nuclear mat- 
ter and negatively charged quark matter with the same 
fj, e coexist at any given radius, with fi e changing with 
radius [19]. If the quark matter is in the CFL phase, 
there is another possibility. At the jl (i.e. the radius) at 
which SIcfl crosses ^nuclear, an interface between bulk 
nuclear matter with N e ^ and bulk CFL quark mat- 
ter with N e = may be stable, as long as the fi e at 
the interface satisfies (12). The CFL insulator cannot 
admit electrons while remaining neutral, even when in 
equilibrium with nuclear matter with large /i e . The de- 
scription of the interface is more complicated than that 
of the bulk phases. Boundary layers within which local 
electric neutrality is not maintained and across which an 
electrostatic potential gradient develops are required, as 
at an ordinary metal-insulator boundary. 
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